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Abstract 

The paper addresses a novel model of metamaterial structure. A 
system of spinners has been embedded into a two-dimensional periodic 
lattice system. The equations of motion of spinners are used to derive 
the expression for the chiral term in the equations describing the dy- 
namics of the lattice. Dispersion of elastic waves is shown to possess 
innovative filtering and polarization properties induced by the vortex- 
type nature of the structured media. The related homogenised effective 
behavior is obtained analytically and it has been implemented to build 
a shielding cloak around an obstacle. Analytical work is accompanied 
by numerical illustrations. 

Keywords: Elastic lattice; chiral material; cloaking. M. Brun et al. Vortex- 
type elastic structure 



1 Introduction 



Dynamic response of vector elastic lattice systems, interacting with waves, 
has been the subject of classical investigations (see, for example, Kunin 1982 
and 1983, Slepyan, 1981), in addition to the more recent publications of 
Slepyan (2002), Slepyan & Ayzenverg-Stepanenko (2002) and Colquitt et al. 
(2011). We note substantial differences between scalar problems involving 
vibration of systems of harmonic springs and vector problems of elasticity, 
referring to elastic rods and beams, which connect a system of finite solids or 
point masses. The notion of shear stress and shear strain is absent in scalar 
lattice systems and there is no such analogue in models of electromagnetism. 
The misconception of predictability of dynamic properties of vector elastic 
lattices is sometimes based on an intuitive extrapolation of results available 
for the scalar systems or for problems of electromagnetism where there is 
only one wave speed. A thorough analysis is the only correct way forward 
and, as illustrated in Colquitt et al. (2011), some results for micro-polar 
elastic lattice structures do not follow from the simpler physical models or 
intuitive assumptions. 

Phononic band gap structures have been analysed in recent years for 
periodic arrays of voids or finite inclusions in the continuum matrix by 
Movchan, Nicorovici & McPhedran (1997), Poulton et al. (2000), Zahpaev 
et al. (2002), Platts et al. (2002) and Lin & Huang (2011). Platts et al. 
(2003)a and Platts et al. (2003)b studied transmission problems for arrays 
of elastic structured stacks including a comparative analysis of the filtering 
properties of elastic waves in doubly periodic media and the transmission 
properties for the corresponding singly periodic stack structure. In particu- 
lar, a review of the Rayleigh Multiple Methods was given by McPedran et al. 
(2001). Here, emphasis is on the control of stop bands in the low- frequency 
region. Furthermore, a highly cited paper (Milton, Briane & Willis, 2006) 
addresses a specially designed elastic coating, with anisotropic inertia prop- 
erties, which may act as an "invisibility cloak", i.e., it routes a wave around 
an obstacle. Brun, Guenneau & Movchan (2008) show an alternative in the 
design of an invisibility cloak which incorporates a micro-polar composite. 

It is also feasible to think of a coating around an elastic inclusion, which 
would enhance (or reduce) the shadow region created as a result of interac- 
tion between the inclusion and an incident wave. Our special interest is in 
the design of a chiral coating, which could be modelled as a homogenised 
micro-structure, characterised by a vorticity constant a. The numerical sim- 
ulation in Fig. [l] shows as interaction of a radially symmetric shear wave, 
created by a time-harmonic point moment, with an elastic inclusion; part 
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Figure 1: Comparison of the cases (a) (a = 0.0) of the elastic inclusion 
without coating and (b) the inclusion with a chiral coating {a = 30.0), 
placed in an elastic medium loaded by a time-harmonic point moment. The 
chiral coating of large a increases the shaded region created by the inclusion, 
and suppresses the displacement. 



(a) of the diagram presents the displacement amplitude in the plane with 
a disk-shaped inclusion without a coating; part (b) shows the high level of 
suppression of the scattered fields when a coated inclusion (with a chiral 
coating) is introduced in the homogeneous medium. The graphs in the bot- 
tom row of the diagram show the amplitude of the displacement measured 
along the main diagonal of the square domain connecting the upper left cor- 
ner to the bottom right corner of the domain. The displacement behind the 
inclusion has been suppressed to a large extent. Needless to say, the appli- 
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cations of this type of design are limitless, with examples including design of 
earthquake resistant structures in addition to drilling systems in geophysics. 

In the present paper, we address propagation of elastic waves, their dis- 
persion properties, and interaction with defects in a two-dimensional struc- 
tured medium containing vortices. Vortices are created by spinning masses 
embedded into the lattice system. The analysis of gyroscopic motion of 
an individual mass is incorporated into the system of conservation of lin- 
ear and angular momenta within the lattice systems. Furthermore, elastic 
Bloch-Floquet waves are considered in such a system. It is shown that the 
presence of spinning masses "stiffens" the overall system, and it also leads 
to a special design of chiral media, which possess stop band properties in 
addition to shielding (re-routing) with respect to elastic waves. In turn, 
we consider several types of homogenised chiral systems to design coatings 
around finite defects in such a way that a dynamic signature of the coated 
defect may be enhanced or suppressed compared to the case when the coat- 
ings are absent. 

2 A vortex-type lattice system 

A general algorithm for analysing the spectral properties of a periodic lattice 
was established in Martinsson & Movchan (2003), and lattice systems with 
built-in dynamic rotational interactions were analysed in Colquitt et al. 
(2011). The new model proposed in the present paper incorporates the 
vortex-type of action at the junction regions of the lattice system. We 
also derive equations containing the vortex terms and explain their physical 
nature. 

2.1 Governing equations 

We consider a triangular periodic lattice with the elementary cell shown in 
Fig. [2j which is defined through the basis vectors: 

t(i) = (2/,0)^ and t(2) = (//2,/V3/2)^. 

In this case, we have two junctions per elementary cell of the periodic struc- 
ture, and it is assumed that the point masses placed at these junctions are 
mi and m2. The stiffnesses of elastic truss-like links, each of length /, are 
assumed to be equal to c, and the mass density of the elastic links is assumed 
to be negligibly small. 
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Figure 2: A triangular lattice linked to the system of spinners is shown in 
part (a). Part (b) shows the geometry of the triangular bi-atomic lattice. 
Different masses at junction points are shown as large and small discs. Part 
(c) shows a spinner and the notations for the coordinate systems; '^/^, and 
9 are the angles of spin, precession and nutation, respectively. 

We also use the three unit vectors = 1,2,3, characterising the 

directions of trusses within the elementary cell: 

a« = (l,Of , a(2) = (-l/2,^/3/2f , a(3) = (-1/2, -^/3/2f . 

The multi-index n = (ni, 77.2) with integer entries is used to characterise the 
position of the elementary cell within the periodic structure, so that the 
position vector of the mass is defined as 

x("''^) =x(0''^)+nit(i)+n2t(2), 

and its displacement time-harmonic, of radian frequency uj. For 

convenience, we also use the multi-indices ei = (1, 0) and 62 = (0, 1). 

The design of the lattice system considered here includes gyroscopic spin- 
ners attached to every junction. The axis of each spinner is perpendicular 
to the plane of the two-dimensional lattice, and its angular velocity can be 
chosen accordingly by comparison with the radian frequency u of the main 
lattice system. A small in-plane displacement of a junction leads to a change 
of orientation of the spinner axis and hence generates a moment creating a 
"vortex-type" effect. 
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The equations of motion for the masses mi and m2 within the elementary 
ceh n are 



,2 

C 

+a(2) . (u(n-ei+e2,2) _ ^{n,l)^^{2) ^ ^_^(2)) . (^(n-e2,2) _ u(n,l))(_a(2)) 

+a(3) . (u("-«2,i) _ u(»'i))a(3) + (-a(3)) . (u("+«2,i) _ u("'i))(-a(3)) 

+^Ru(--i), (1) 

c 

and 

_!!!H^u(n,2) ^ ^(1) . („(n+e„l) _u(n,2))^(l) + . „(„,2) ) ) 

c 

■ (u("+^2'l) - u("'2))a(2) + (-a(2)) . (u(n+ei-e2,l) _ u(n,2))(_^(2)) 
+a(3) . (u(n-e2,2) _ u(n,2))^(3) + ^_^(3)) . ^^{n+e„2) _ ^(",2) ) ^_^(3) ) 

+ ^Kn(-''\ (2) 

c 

where R is the rotation matrix 
and ai, q;2 are the spinner constants. 



2.2 Evaluation of the spinner constants 

It is assumed that each junction within the lattice is connected to a spinner, 
whose axis is perpendicular to the plane of the lattice, as shown in Fig. 
[2| By considering one of the junctions, together with the spinner attached, 
we evaluate the spinner constants aj^j = 1,2, introduced in the previous 
section. 

Let ^, (/) and be the angles of spin, precession and nutation with respect 
to the vertical axis Oz, respectively (see, for example, Goldstein et al., 2000, 
page 210). An individual spinner and the angles ^, and are shown in 
Fig. Then, assuming that the gravity force is absent, the equations of 
motion of the spinner, for the case of the constant spin rate = 0), are 

= /o {0 - (j)^ sin cos 0) + sin 0{(j) cos 6> + ^) , 
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My = sin e [I^ {(/) sin 6 + 2(1)0 cos 0) - 10{(j) cos 6> + ^) j , 

M^ = I{4) cos 0-^0 sin (9) , (3) 

where Mx, My, are the moments about the x,y and z axes, and /q = 
Ixx — lyyi I — Izz circ the moments of inertia. 

We consider a smah amphtude time harmonic motion (due to the motion 
of the lattice attached to the spinner), resulting in the nutation angle 

0{t) = ee^^^ |e| < 1, 

and constant spin rate ip = Q. 

Assuming that there are no imposed moments about the x and y axes 
(i.e. Mx = My = 0), and the precession rate is constant (0 = 0), we deduce 
the connection between the spin rate ft and the radian fequency uj of the 
imposed nutation. Namely, the precession rate is given as 

■ (4) 



" 2Io - / 

According to the equations of motions we have 

(/ - /o)02 + - Iqu^ = 0. (5) 
The compatibility condition for equations Q and ^ has the form 

Direct substitution into ([s]) gives the induced moment about the z— axis 

M, = ^Iuj00. 

Let h be the characteristic length of the spinner, as shown in Fig. Then 
the magnitude of the in-plane displacement in the lattice junction is 

U = 0h. 

Taking into account that the moment about the z— axis is Mz = Fh0, we 
deduce that the "rotational force" F is given by 

Following the above derivation, the spinner constants ai and 0:2 in equations 
^ and Q are 
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2.3 Elastic Bloch waves 



Consider the case when mi = 1712 and o^i = 0^2, i-e. the mono-atomic lat- 
tice. Then the equations ([T]), ^ are equivalent, subject to a translation 
of an elementary cell of the periodic system, which includes only a single 
junction. This configuration is a good example for observing the influence 
of the vortex-type interaction within the lattice on the dispersion properties 
of the elastic waves. The amplitudes Ui, U2 of the time-harmonic displace- 
ment, of radian frequency cj, within the periodic vortex-type lattice satisfy 
the system of equations ([T]) together with the Bloch Floquet conditions set 
within the elementary cell of the periodic system, as follows: 



u(x + nit(^V2 + ^2t 



u(x) exp(zk • Tn) 



(6) 



where k is the Bloch vector (fci, k2), and the matrix T has the form 

I 1/2 

1V3/2 



These waves are dispersive, and the corresponding dispersion equation 



IS 



det 



C(k) -w2(M 



0, 



(7) 



where the mass matrix M = diag[m, m] and SI represents the chiral term 
associated with the presence of spinners attached to the junction points 
of the lattice system. Assuming that each spinner is characterised by the 
"spinner constant" a, we represent the matrix SI as 



—ia 
ia 



(8) 



Additionally, in equation Q, C(k) is the stiffness matrix for the monatomic 
lattice of masses m, with spring connectors of stiffness c. It is given by 



C(k) 



/„ „ , , (C0S$ + C0S^') \/3(cOS * - COS $) \ 

'3 — 2 COS kil — ^ ' 



V^(C0S * - COS $) 



3(cos$ + COS*) 



where 



$ = — + —k2l and ^' 



kil v^, , 



We note that the dispersion equation is bi-quadratic with respect to cj, 
and it has the form 



8 



^^(m^ - a^) - u?m trC + detC = 0. 



(9) 



Since c > 0, both trC and detC are positive for all ki and k2 in the 
elementary cell of the reciprocal lattice, except at the origin where they are 
both zero. Hence there are two important regimes. For rn? > o? ^ there are 
two dispersion surfaces. As the factor {m} — o?^ decreases towards a critical 
point where — o? ^ the upper dispersion surface increases without limit 
and the lower surface decreases. At this critical point, the upper surface is in- 
finite and the lower surface has the degenerate value of {m~^ detC/trC)-*^/^. 
For VP? < a^, the lower dispersion surface only remains. Representative 
results will now be given for each of the regimes > (subcritical) and 

< o? (supercritical). 

2.4 Dispersion properties 

Firstly, Fig. |3] gives dispersion surfaces and the corresponding two-dimensio- 
nal diagrams with the slowness contours for the uniform lattice without spin- 
ners (q: = 0). Two conical surfaces are clearly visible in the low frequency 
range. These surfaces are shown to evolve as the spinners are brought into 
the system, now with chiral terms in the equations of motion. 

2.4.1 Monatomic lattice of the vortex- type 

The diagrams for the vortex-type of monatomic lattice, with a small spinner 
constant!^ a = 0.5 and mass m = 1 are given in Figjij It is apparent that the 
lower acoustic surface, dominated by shear waves, does not change by much 
as a result of the vortex interaction. On the contrary, the upper acoustic 
surface, dominated by pressure waves extends further into a higher range of 
frequencies and has a pronounced conical shape for a vortex-type of lattice 
of small a. The latter implies that the homogenisation range for Bloch- 
Floquet waves corresponding to this dispersion surface is also extended into 
a wider range of frequencies, as shown in Fig. |4j 

For the other regime, the dispersion diagram illustrating such a situation 
for the case of a = 2 and m = 1 is given in Fig. [5j Only the single dispersion 
surface is seen for the low frequency range. 

^Normalisation has been used throughout the text, so that c — 1^ m — 1 and the 
distance between the neighbouring masses is equal to unity. All other physical quantities 
have been normalised accordingly, and the physical units are not shown. 
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The influence of the vortex-type of interaction appears to be significant 
for the dispersion properties of Bloch Floquet waves within the entire admis- 
sible frequency range. The illustration is given by Fig. |6] where we assume 
ki — k2 — k and plot c<j as a function of the spinner constant a and the Bloch 
parameter k. In this case, the lower dispersion surface is represented by the 
diagram in Fig. [6^ and is defined for all real values of a whereas the up- 
per dispersion surface disappears in the supercritical regime when \a\>m. 
The diagram characterising the upper dispersion surface and showing uj as 
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(a) 




(b) 

Figure 4: (a) Dispersion surfaces for the monatomic lattice with spinners a 
= 0.5. (b) Contours for the two eigenfrequencies for the monatomic lattice 
with spinners o: = 0.5. 



a function of a and k is presented in Fig. [6)3 for < o; < m where m — 1. It 
is noted that for a fixed value of A:, the frequency uj on the upper dispersion 
surface increases with the increase of the spinner constant a (see Fig[6]3), 
whereas on the lower dispersion surface the frequency uj decreases as the 
spinner constant a increases (see Fig. lofe 



11 



kll 



kol 



Figure 5: Dispersion surface for the monatomic lattice with spinners a = 2. 



2.4.2 The low frequency range 

We examine the behaviour of a;(k) in the low frequency limit for small values 
of |k|, and evaluate the effective group velocity corresponding to a quasi- 
static response of a homogenised elastic solid. The solutions to equation Q 
may be expanded for small values of ki and k2 as follows 



'^l = n/-((*^l')' + (W). a = m. (11) 
8 V m 



and 



/3c /2m + (m2 + 3a2)i/2\ 
For the case when no spinners are present, these formulae reduce to 



- • / ^ {{kiiy + {k2lY) and ujI = V3 c^?. 



These expressions for uj^ and consistent with the special case of 

equation (3.3) in Colquitt et al. (2011). 
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(a) 



(b) 



Figure 6: The radian frequency change at the lower dispersion surface (a) 
and the upper dispersion surface (b) with a along the line ki = k2 in the 
reciprocal space, m — 1^ c — 1 and / = 1. 



In the subcritical regime m? < o? ^ for low frequencies, only one dis- 
persion surface exists representing a shear wave. For a monatomic, har- 
monic scalar lattice (with no spinners), the dispersion equation is given in 
Ayzenberg-Stepanenko & Slepyan (2008) as 



(jj — 



\ 



o kil ^ kil V3k2l 

4 cos^ 4 cos — cos 

2 2 2 



where the subscript ^corresponds to the scalar case of out-of-plane shear. 
The low frequency approximation for this dispersion equation is 



3cg 
2m. 



(12) 



It is interesting to compare equations (10), ( |11| ) and (12) in the regime 
< a^. Because of the presence of only one dispersion surface in this 
regime for the vector lattice of elasticity, then this may be compared with 
an 'equivalent' scalar lattice in which the waves have the same dispersion 
properties as the shear waves in the vector lattice in the low frequency 



regime. Comparing equations (10) and (12) and by the appropriate choice 
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Cs _ c v^l + 3(a/m)2-2 

rris 4m {a/mY — 1 

of combination of parameters Cg, c and m^, m for both lattices together with 
a for the spinners in the vector lattice, the shear wave in the vector lattice 
and the wave in the scalar lattice may be constructed so as to have the same 
dispersive properties. 

2.4.3 Bi-atomic lattice of the vortex-type 



mi =1 7712 =10 c =1 I = \ Qi = q;2 = 



mi — 1 7712 — 10 c— 1 I — 1 OLi — 0.5 a2 — 0.5 




(a) 



(b) 





kil = /C2/ 



kil — 



(c) 



(d) 



Figure 7: Dispersion surfaces for the biatomic lattice (a) with contrasting 
masses and no spinners, and (b) with contrasting masses and the same 
spinner throughout the lattice, (c) and (d) are the cross section of the 
surfaces given in (a) and (b), respectively, for kil — k2 I. 



Now we show the computations for Bloch-Floquet waves for the biatomic 



triangular lattice described in Section 2.1, The dispersion equation is given 
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by relation Q, where 



M = diag[mi, mi, m2, 7112] 

and 



5](a) 



/ -iai ' 
iai 
-ia2 
V ia2 / 

Additionally, the stiffness matrix can be expressed in the form (see Martins- 
son & Movchan 2003) 



C(k)= c 
with the 2x2 matrices 

Cii(k)=C22(k) = 



Cn(k) Ci2(k) 

C2l(k) C22(k) 



3 - I cos $ cos$ 
-^cos$ 3 - |cos$ 



C2i(k) = Ci2(k) = exp[z($ + 



-2 cos($ + - I cos ^ cos 

^cos^^ -|cos^^ 



For such a lattice with no spinners, the dispersion surfaces for Bloch- 
Floquet waves are shown in Fig [7^ (a section of the surfaces for ki — k2 
is added in Fig. [t]^). In contrast, for the vortex- type of lattice the cor- 
responding results are given in Fig (with a section for ki — k2 in Fig. 

These computations show that the vortex-type of interaction within the 
lattice leads to formation of additional total band gaps, as shown in Fig. 
[7)3. It is noted that the change in the behaviour of the lowest two dispersion 
surfaces through the introduction of spinners is small, whereas additional 
band gaps are introduced in the upper two dispersion surfaces. 

The influence of the spinner constant a on the dispersion properties of 
Bloch-Floquet waves is illustrated in Fig. |8] where u is represented as a 
function of a and kl = kil = ^2/. The two lower shear dispersion sur- 
faces, represented in Fig. [8^, are defined for all real values of a, whereas 
the two upper pressure dispersion surfaces, depicted in Fig. [8)3, disappear 
when critical regimes determined by the spinner constant a are reached. In 
particular, in the subcritical regime {a < 1) two pressure waves propagate 
while one disappears in the intercritical regime (1 < a < 10) and no pressure 
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m-i =1, 1712 =10, c =1, / =1, ai — a2 — Q: 



mi =1, 1712 =10, c =1, I =1, ai = q;2 = « 




/ci/ = k2l 



kil — k2l 



(a) 



(b) 



Figure 8: The radian frequency change at the two lower dispersion surfaces 
(a) and the two upper dispersion surfaces (b) with a along the line ki = /c2 
in the reciprocal space, mi — 1, 1712 = 10, c = 1 and / = 1. 

waves propagate in the supercritical regime (a > 10). Also, for the lower 
dispersion surfaces, at sufficiently large a the frequency u decreases as the 
spinner constant a increases (see Fig. Isk). 



3 Homogenisation approximation. Dynamic shield- 
ing 



Here we discuss the dynamic response and shielding properties of a ho- 
mogenised elastic solid, whose equations of motion contain the vortex-type 
term. 

As illustrated in the previous section the long- wave approximation for the 
triangular lattice, incorporating the periodic system of spinners, corresponds 
to an isotropic elastic material with the Lame constants X — ji^ mass density 
p and an additional term representing the chiral properties of the medium. 
The amplitude vector U of the time-harmonic displacement satisfies the 
equations of motion as follows 



where /i is the shear modulus of the homogenised solid, SI is the vorticity 
matrix defined in ([s]). In the computations presented here, we normalise 
the mass density and the shear modulus of the chiral medium, so that p — 



AU + 2VV • U ) + cj^EU + pcj^U + F = 0, 
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1.0, /i = 1.0, and distinguish between (a) the subcritical case of \a\ < 1 and 
(b) the supercritical case of\a\ > 1, where a is the spinner constant from ([s]). 
Such a classification is hnked to the dispersion equation Q. In the regime 
(a), there are two conical surfaces for small values of u in the dispersion 
diagram (see Fig. [4]). On the other hand, in the regime (b) there is just one 
acoustic dispersion surface, which has a conical shape for small values of u. 

The chiral medium is used in the coating for an elastic inclusion, which 
interacts with an incident wave created by point sources of different types. 



(a) (b) 




(c) (d) 

Figure 9: Displacement amplitude produced by the unit force applied in the 
xi— direction in the homogenised chiral medium: (a) a = 0.0 (no spinners), 
(b) a 0.6, (c) a 1.0, (d) a 1.4. 
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3.1 Point force in the chiral solid 

Firstly, we show the results of the computational experiment produced for 
a homogenised chiral solid loaded by a time-harmonic horizontal point force 
of the radian frequency = 50. 

Fig. [9] shows the change in the wave pattern for different values of the 
spinner constant a. In particular, Fig. [9^ has the result of the computation 
for the non-chiral medium {a — 0.0), whereas Fig. ^jp includes the modified 
plot corresponding to a subcritical value of the spinner constant {a — 0.6) 
demonstrating a vortex around the point force. Domination of pressure 
waves, with a relatively large wave-length is clearly visible in the horizontal 
direction in Fig. [9^ compared to Fig. [9)3. 

Fig. ^ shows the critical case, with the spinner constant a = 1.0, and 
the presence of spinners leads to a nearly isotropic map; the preferential 
directions associated with the point force are less pronounced. 

Furthermore, Fig. [9]l includes the displacement amplitude for another 
supercritical case with the larger value a = 1.4, and the dominance of shear 
waves, with the relatively small wave length is clearly visible. 

The computations suggest that the directional preference is substantially 
reduced in the presence of 'spinners' in the chiral medium in the supercritical 
regime of \a\ > 1. 




Figure 10: Coated inclusion. The coating consists of the chiral material with 
the spinner constants a or —a, as illustrated on the diagram. 
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3.2 An elastic inclusion with a chiral coating 



Here we consider an elastic plane containing a coated inclusion. The two- 
phase coating of the inclusion is assumed to have the chiral term in the 
governing equations, whereas the ambient medium, around the coating and 
the inclusion itself are assumed to be non-chiral. The coating consists of 
two semi-rings, with the spinner constants being of the same magnitude and 
opposite sign, as shown in Fig. 
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The coated inclusion is chosen to be 
symmetric about the straight line passing through the centre of the inclusion 
and the point where the concentrated force is applied. 




(a) (b) (c) 

Figure 11: Displacement amplitude in the plane containing a coated inclu- 



sion (as shown in Fig. 10). The loading is provided by the horizontal point 
force. The inclusion is placed in the region dominated by the shear waves. 
Three computations correspond to different values of the spinner constant 
a: (a) a = 0.0; (b) a = 0.9; (c) a = 2.0. 



An incident wave is generated by a point source placed at a finite distance 
from the inclusion. For the computational examples, the source term is 
provided by a point force or concentrated moment with a time harmonic 
amplitude of radian frequency cj = 50. In the computation, the elastic 
inclusion has the normalised Lame constants A = 23, ^ = 12 and the same 
mass density as in the ambient medium. The coating is assumed to have 
the same elastic moduli and the same mass density as in the ambient elastic 
matrix. 

For the isotropic medium with an inclusion without a coating {a = 0.0) 
the results of the computation are shown in Fig. [TT^. Comparing with Fig. 
[9] one can clearly see the shaded region associated with the perturbation 
field produced by the elastic inclusion. 

The next set of computations in Figs. [TT] 3 and [TT] ^ corresponds to the 
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(a) 



(b) 



Figure 12: Comparison of the cases (a) {a = 0.0) of the elastic inclusion 
without coating and (b) the inclusion with the chiral coating (a = 1.5), 
placed in the elastic medium loaded by the time-harmonic point force 
pointed towards the centre of the inclusion. The chiral coating reduces 
the shaded region created by the inclusion. 



chiral two-phase coating placed around the inclusion and can be interpreted 
as a "cloak" guiding shear waves around the inclusion. It is noted that, 
in these two diagrams where the point force is horizontal, the inclusion is 
placed in the region dominated by shear waves. It is also seen that the effect 
of cloaking is most pronounced in the supercritical regime when the spinner 
constant of the inclusion is chosen so that \a\ > 1. 



One of the important features of the computations shown in Fig. [TT 
is the directional preference of the ambient field produced by the point 
force. In particular, the above diagrams correspond to the case when the 
ambient field in the immediate neighbourhood of the inclusion is dominated 
by shear waves. A change of the orientation of the force with respect to the 
position vector of the inclusion results in a different pattern of the ambient 
field surrounding the inclusion. Fig. [12] correspond to the situation when 
the force is pointed towards (or opposite) the centre of the elastic coated 
inclusion. As before we show the results of the computations for the inclusion 
without the coating [a = 0) in Fig. 12 a., and supercritical chiral coating 
(a = 1.5) in Fig. [T2| 3. The direct comparison of these diagrams shows that 
the chiral coating reduces the shadow region associated with the inclusion. 
Furthermore we observe that the ripples behind the inclusion have a size 
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(a) (b) 

Figure 13: Comparison of the cases (a) homogenous medium with no inclu- 
sion and (b) the inclusion with the chiral coating {a = 2.0), placed in the 
elastic medium loaded by the time-harmonic point moment. 

corresponding to that of the wave length of shear waves, which suggests 
that the chiral coating is acting as a shear polarising cloak. 

Finally, we consider the configuration when the ambient elastic field 
has no directional preference. The loading is provided by a time-harmonic 
concentrated moment of the radian frequency u = 50. The computation for 
the homogeneous medium, without an inclusion is presented in Fig. [13^. 
We also note the earlier Fig. [l^ where the computation is presented for an 
uncoated inclusion. The shadow region can be reduced via introduction of 
the chiral coating. The supercritical case of o: = 2.0 is included in Fig. [13) 3, 
which shows a new diffraction pattern, with conical regions where waves have 
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enlarged amplitude. Two diagrams in the bottom row of Fig. 13 give the 
amplitude of the total displacement along the main diagonal of the square 
computational region. They show that along this line the displacement 
amplitude behind the inclusion with the chiral coating of a = 2.0 is similar 
to the displacement amplitude in the homogeneous medium without any 
inclusion. 



4 Concluding remarks 

This paper has placed together two concepts: the model of a discrete lattice 
with spinners and the homogenisation approximation by a continuum chiral 
medium. The vorticity constants in the governing equations describing the 
chiral media are evaluated explicitly and, furthermore, dynamic response of 
both discrete and continuous systems have been analysed in detail. 

Novel properties of Bloch waves in discrete vortex-type elastic lattices 
have been identified, with analytical findings being complemented by the 
numerical illustrations of dispersion surfaces. 

In the homogenisation approximation, the chiral material has been used 
to design a composite cloak around an inclusion, so that the incident waves 
are guided around the inclusion. The result of such an interaction can lead 
to a substantial reduction of the shadowed region behind the inclusion, and 
the observation that the coating can be interpreted as a polarising cloak. 

Furthermore, it was demonstrated that the parameters of the chiral coat- 
ing can also be tuned so that the shadowed region is enhanced and the am- 
plitude of the displacement behind the inclusion becomes negligibly small. 

We envisage a range of applications of the proposed model in problems of 
geophysics and structural design of cloaks shielding defects, which interact 
with elastic waves. 
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